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AN EXTENSION PROBLEM FOR THE FRACTIONAL DERIVATIVE 

DEFINED BY MARCHAUD 

CLAUDIA BUCUR AND FAUSTO FERRARI 


Abstract. We prove that the (nonlocal) Marchaud fractional derivative in R can be 
obtained from a parabolic extension problem with an extra (positive) variable, as the 
operator that maps the heat conduction equation to the Neumann condition. Some 
properties of the fractional derivative are deduced from those of the local operator. In 
particular we prove a Harnack principle for Marchaud-stationary functions. 
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1. Introduction 


In literature there are several definitions of fractional derivatives (see, for instance, the 
monographs [20], [19] and [4] for an historical introduction). In particular, we are interested 
in the notion given by Marchaud, see [16], who introduced two types of fractional derivatives. 
For a fixed s € (0,1), the left and the right Marchaud fractional derivative of order s, see 
[20], formulas 5.57 and 5.58, are respectively defined as follows: 


Di/M 


_ r m-f(t tt ) 

r(i - s) Jo T l +‘ 


( 1 . 1 ) 


These fractional derivatives are well defined when / is a bounded, locally Holder continuous 
function in 1R. In particular, we may assume that / £ C^(M), for s < 7 < 1 and / £ L°°(R) 
(see the Appendix for further details), even though these hypotheses can be weakened. In 


Key words and phrases. Fractional derivative, Marchaud derivative, extension operator, Harnack in¬ 
equality. MSC: 26A33, 35K10, 35K65. 
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addition, we just recall here that the Marchaud derivative can be defined for s £ (0, n) and 
n £ N, as 




{4 r f [s] (t)-fW(t T T) 

r(i -{s}) J 0 n+w 


where [s] and {s} denote, respectively, the integer and the fractional part of s. Our work 
focuses on the case n = 1 and, in the first part of the paper, on the left fractional derivative, 
that we can write by a change of variable, neglecting the constant and omitting for simplicity 
the subscript symbol +, as: 


d s m 


JO T S+ 



f(t) - Hr) 

(t — t) s+1 


dr. 


( 1 . 2 ) 


A short remark on the right counterpart of the Marchaud fractional derivative is given 
in Section 5. Moreover, we point out that we argue considering (1.2) as the definition of 
our fractional derivative without taking care of what happens when s —> 0 + or s —> 1~. 
Nevertheless, in the Appendix, we briefly discuss these interesting cases with respect to the 
definition given in (1.1). 

The purpose of the present work is to introduce an extension operator for the fractional 
derivative introduced in (1.2). Indeed, the operator D s naturally arises when dealing with 
a weighted parabolic differential equation (the heat conduction problem) on the positive 
half-plane, with a positive space variable and for all times, namely for (x,t) £ [0,oo) x R. 

In order to construct this extension operator, we exploit the idea recently revisited in 
[5]. In that paper, the fractional Laplacian was characterized via an extension procedure, 
by means of a weighted second order elliptic local operator. 

Considering the function p of one variable, formally representing the time variable, our 
approach relies on constructing a weighted parabolic local operator by adding an extra 
variable, say the space variable, on the positive half-line, and working on the extended 
plane [0, oo) x M. 

The heuristic argument can be described as follows. Let p : R —> R be a given function, 
sufficiently smooth. Let U be a solution of the problem 


dU _ d 2 U 
~dt ~ 0a?’ 
U(0,t) = <p{t), 


(, x , t) £ (0, oo) x R 
t £ R. 


(1.3) 


Let us point out that this is not the usual Cauchy problem for the heat operator, but a 
heat conduction problem. 

It is known that, without extra assumptions, we can not expect to have a unique solution 
of the problem (1.3), see [21]. Nevertheless, if we denote by T) / 2 the operator that associates 


dU 

to p the partial derivative , whenever U is sufficiently regular, we have that 

ox 


dp 

dt 


T1/2T1/2P 

That is Tj /2 acts like an half derivative, indeed 








AN EXTENSION PROBLEM FOR THE FRACTIONAL DERIVATIVE DEFINED BY MARCHAUD 3 


The solution of problem (1.3) under the reasonable assumptions that p is bounded and 
Holder continuous, is explicitly known (check e.g. [21]) to be 


U(x,t) = cx / e 4 u- t ) (t — t) 2 p(t)cIt 


L 


—oo 
oo 


= cx e 4t t 2 p(t — t) dr, 


x 2 

where the last line is obtained with a change of variable. Using t = — and the integral 

4r 

definition of the Gamma function (see formula 6.1.1 in [1]) we have that 


dr = 2 


e *t 2 dt = 2T 



Hence, 


U(x,t) - U(0,t) 


x 


= c / e 4t r 2 (p(t — r) — p(t)) dr , 

J o 


choosing c that takes into account the right normalization. This yields, by passing to the 
limit, that 

nOO 


U(x,t) — U(0,t) 

— lim —- = c 

x— > 0 + X 


p{t) - p{t - t) 


dr. 


o 


r 2 


(1.4) 


Hence, with the right choice of the constant, we get exactly T) l t 2 ip i.e. the Marchaud 
derivative of order 1/2 of p. 

Now we are in position to state our main result. 


Theorem 1.1. Let s £ (0, 1) and p: M — > R be a bounded, locally C 7 function for s < 7 < 1. 
Let U: [0, 00 ) xIaR be a solution of the problem 


dU(x,t) 1 — 2 sdU{x,t) 

dt x dx 

U(0,t) = p(t), 


lim U(x, t) = 0. 

x—>oo 


d 2 U(x, t) 
dx 1 


( x , t) £ (0, 00 ) x R 
t £ R 


(1.5) 


Then U defines the extension operator for p, such that 

D s p(t) = - \im + c s x~ 2s (U(x,t) - p(t)), ( 1 . 6 ) 

where 

c s = 4 s T(s). 


We notice that one can write 

f)TI 

D s p(t) = - lim CsX 1 ^ 28 — (x,t), (1.7) 

x->-o+ dx 

in analogy with formula (3.1) in [5]. 

Remark 1.2. The extension defined in (1.7) satisfies, as one would expect, up to constants: 

D^D s p{t)=p'(t). 
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Indeed, 



( 1 . 8 ) 


An interesting application that follows from this extension procedure is a Harnack in¬ 
equality for Marchaud-stationary functions in an interval JC1, namely for functions that 
satisfy D s y; = 0 in J. This result is not trivial, since the fractional stationary functions on 
an interval of R determine a nontrivial set of functions, see e.g. [3]. 

Theorem 1.3. Let s € (0, 1). There exists a positive constant 7 such that, ifD s ip = 0 in 
JCI and ip > 0 in R, then 


(1.9) 


sup tp < 7 inf ip 
[t 0 -f 5,t 0 -j5] [to+f <5,to+<5] 

for every to £ R and for every 8 > 0 such that [to — 5, to + 5] C J. 


This result can be deduced from the Harnack inequality proved in [G] for some weighted 


parabolic operators. In particular, the constant 7 used in the previous Theorem 1.3 is the 
same that appears in the parabolic Harnack case, see [6]. In addition, we remark that 
the inequality (1.9) is not the usual Harnack inequality for elliptic operators, where the 
comparison between the supremum and the infimum is done on the same set, e.g. the same 
metric ball. This Harnack inequality for the Marchaud-stationary functions inherits the 
behavior of its parabolic extension. 


2. The weighted parabolic problem 

In this section we find a solution of the system in (1.5). At first, we introduce a particular 
kernel, that acts as the Poisson kernel. We then look for a particular solution of the system 


by means of the Laplace transform, and in this way we show how the solution arises. Finally, 
by a straightforward check, it yields that indeed the indicated solution satisfies the problem 


(1.5). 


2.1. Properties of the kernel In this section we introduce and study the properties 
of a kernel, that acts as the Poisson kernel for the problem (1.5). 

We define for every ieK, 



( 2 . 1 ) 


Also, let 



( 2 . 2 ) 


and notice that 



(2.3) 
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Indeed, we have by changing coordinates r = —^ 

ar 

L*- {x ’ t)dt = mr)i 


x 2s e~^ T t~ s ~ 1 


dt 


4sr (s) Jo 

= f ip s (t)dt. 

Jr 

The kernel satisfies also the following property: 


e 4t t 


dr 


f ^ s (x,t) dt = 1. 

JR 


Indeed by recalling (2.3) and performing the change of variables t = —, we get 

4 r 


/ 

Jr 


ips (t) dr = 


r °° _L- _ a _i , 1 

e 4 t T 


dr = 


e~ T T s ~ l 


dr. 


4s r(s) Jo r(s) 

From the integral definition of the Gamma function (see formula 6.1.1 in [1]) 


(2.4) 


(2.5) 


e t t s 1 dt , 


it follows that 


F(s) = f 

do 

[ *1 ’ s {x,t ) dt = 1. 

JR 


Taking the Laplace transform of the kernel \H S , we have the following result involving the 
modified Bessel function of the second kind K s , see [15] and [1], §9.6. 

Lemma 2.1. The Laplace transform of the function ip s £ L 1 (M.) is 


1 


C(lp s ){Lo) = ———— W = K s(y/uj) for SfttJ > 0. 


( 2 . 6 ) 


2 s_1 r(s) 

Moreover, the Laplace transform with respect to the variable t of the kernel Ts £ L 1 (R, dt) 
is 


£(4' s )( x,u) = 1 x s u}*K s (xy/w) for > 0. 

2 s X 1 (s) 


(2.7) 


Proof. If one proves claim (2.6), the second result (2.7) follows after a change of variables. 
We have that 

1 


£(^g)(a;,w) = 


4 s T(s) 


s 

Jo 


J2s ^rr 4- — S— 1 


Taking t = x 2 r (and recalling that x > 0), we obtain 


£(4' s )(x,w) = 


S -^T - s - 1 e~ u( ~ x 2 t Ut 


4 s r(s) J o 

= £(i/> s )(x 2 w). 

For > 0 and ui £ C with > 0, as stated in formula 5.34 in [18], we have that 


C 


(tP 1 e °) = 2 ^ 2 K 7 ^2(aw)^ . 
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Taking 7 = —s, recalling that K s = K_ s , we obtain 

/> OO 

/ e-^r- 8 - 1 e-“ T dr = 2 s+1 w^K s (V^). 


Hence, multiplying by 4 s T(s) we get: 

1 1 


4 s T(s) 


e -TF r _s_ 1 e _ “ T 


1 


dT= 


2 s - 1 r( s ) 

and therefore ( 2 . 6 ). □ 

2.2. Existence of the solution. We prove in this section the following existence theorem: 
Theorem 2.2. There exists a continuous solution of the problem (1.5) given by 

U(x,t) = ip(t) := f 8 {x,T)<p(t-T) dr. (2.8) 

Jr 

More precisely (inserting the definition (2.1)) we have that 

1 


U(x, t ) = 


„2 s 


4 s T(s)‘ 


e 4t t s ip(t — t) dr. 


(2.9) 


Before proving this theorem, we recall a useful result (see [9], Proposition 4.1) involving 
the modified Bessel function of the second kind. 


Proposition 2.3. If — 00 < a <1, the boundary value problem 

x a y"(x) = y(x) in (0, 00 ) x R 
y(0) = 1 
lim y(x) = 0 . 

x—*oo 

has a solution ip € C 2 ~ a ([0, 00 )) of the form 


ip(x) = c k x 2 Kj_ 

V k 


where 


is a positive constant and 


Ck = 


k := 


2 1 2k k 2k 

r wr 

2 — a 


( 2 . 10 ) 


( 2 . 11 ) 


We continue by showing how the solution of problem (1.5) arises, by means of the Laplace 
transform (see [7] for details on this integral transform). 


Proof of Theorem 2.2. First, we look for a possible candidate of a solution of the problem 
(1.5). In particular, we put ourselves in the simplified situation that U has a sub-exponential 
growth in t and that the function <p is zero on the negative semi-axis (— 00 ,0]. Under this 
additional hypothesis, we take the Laplace transform in t of the system (1.5). Since the 
derivative of the Laplace transform acts as 

£(/')M = u£f(u), 
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we get that 


m 

in 


rTT , . 1-2 sdCU. , d 2 £U. , . /n , 

uj£U(x,uj) = —-- + 2 (x,o>), m (0,oo) x 

£1/(0, t) = 
lim £U{x,uS) = 0, 

tc—>-oo 

We define 

/0*0 := £U(x,uj), 

then / must be solution of the system 

w/(z) = -— —f'(x) + f"{x), in (0, oo) 


/( 0 ) = Cip{u) 
lim f(x) = 0 . 

x—>oo 


( 2 . 12 ) 


(2.13) 


We consider the problem (2.10) for a = (notice that for s £ (0,1), we indeed have 
that a € (—oo, 1 )) 

x~‘~ y"(x) = y{x) in ( 0 , oo) x K 
2 /( 0 ) = 1 
lim y(x) = 0 . 

v x—too 

We observe that by taking 

f(x) = £<p(u)y 

one obtains that / satisfies problem (2.13). Indeed, we have that 

(“• (sf) - 


, 2 s 


(2.14) 


and 




1 


uj~ 2s (2s) 4s ~ 2 x 2 ~ 4s 


£ip(w 
1 - 2 s 


+ /'(*)-^T^) 48 - 2 ^- 2 ^ 1 - 4 *. 

£tp(u)) 


(A' 2 ' 

\ 2 s 


x \ 


2s 


y " \ uS [^) ) =y\ u > t 


(s) 


2s 


Then, since by (2.10) we have 

;) 

this implies that 

/(x) = a ;- 1 (f”(x) + (1 - 2 s)x _ 1 /'(x)) , 
which is the first equation in (2.13). The fact that /(0) = £(p(w) easily follows. 

A solution of problem (2.10) is given explicitly in Proposition 2.3. For a = 28-1 , we have 


k = tt and 

2s 


2 1 _ 8 ( 2 s ) 8 i 


Inserting this into (2.14), we obtain that 


2 / 0*0 = / x^K, 

r(s) 

that 

2 1-s 


( 2 sx 2s ^ . 


(2.15) 


/(x) = £tp(uj)——U}2X S K s (xy/uj). 

r(s) 
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Hence, by (2.12), we have that 

2l-S 

CU(x,Uj) = £ip(uj) . UJ^X S K s (xy/uj). 

r(s) 

By taking the inverse Laplace transform, recalling that the pointwise product is taken into 
the convolution product, we obtain that 

U(x,t) = -ip* £ _1 (u%x s K s (xy / w)) (t). (2.16) 

r( s ) 

From (2.7) we have that 

£(^ s )(a;,w) = (xy/w). 

Hence in (2.16) we obtain the following Laplace convolution 

U(x,t) = <£* ( 2 -l?) 

J 0 

We recall now that we took the function ip to be vanishing for t £ (— 00 ,0). Hence, it 
is reasonable to suppose that the above formula holds true also for a function that is not a 
signal. Hence, now we assume that p is defined on the entire axis R, without the assumption 
that ip vanishes in (— 00 , 0). We claim that p* 'F s still defines a solution of the problem (1.5) 
as in (2.8). Indeed, we show now that the function U defined in (2.8) solves the system 
(1.5). 

In definition (2.8), we call 


and notice that 


A 


X,T • 


e 4t r s 1 , if t > 0 
0 if r < 0 


dA x 


dx 


X 

—— A x r , if r > 0 
2 T ’ 

0 if t < 0. 


Let 


V(x,t) := 4: S T(s)U(x,t) = x 2s f A x>T p(t — t) dr. 

Jr 


Taking the derivative with respect to x of V ( x , t) we have that 

2s+l /* A 

I -ft-cr. 


dV 

dx 


I" T ^s+l r A 

(x, t) = 2sx 2s ~ 1 / A x<T p{t — t) dr — — / —^p(t - r) dr , 

Jr 2 j r t 


and 


dx 2 


d2V (x,t) = 2s(2s - l)x 2s ~ 2 [ A XjT p(t — r)dr - f - r) dr 

Jr 2 Jr t 

™2s+2 r A 

+ — ^f<P(t-T)dT. 

4 JR T 


Then, by changing variables, we write 


V(x,t) = x 2s 


A x , t -T^p{x) dr, 
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and taking the derivative with respect to t we have 


dV 

~dt 


0 x,t) = x 2s [ 
JR 


) A-X,t — T 

4 (t — r) : 


Mr) + 


<p{t) 


dr. 


We change back variables to obtain 

(x, t) = -x 2s+2 j <p(t — t) dr — (s + l)x 2s J — ^v?( t - t) dr. 

By substituting these computations, we obtain that indeed V, hence U by the definition of 
V, satisfies the equation 


dU{x,t) 1 — 2 sdU(x,t) d 2 U(x,t) 
dt x dx dx 2 

Moreover, using for x large enough the bound 

x 2 1 

x 2s e~— < Me -17 , 

thanks to the Dominated Convergence Theorem and the limit 

x 2 

lim x 2s e~~ = 0, 

x—too 


it yields that 


lim U(x, t) = 0. 

x—too 


Furthermore, in (2.8) by changing variables f = —~ (but still using r as the variable of 

x z 

integration), we have that 

1 f°° i 

U{x,t) = 4ar , ' J e~^T~ s ~ 1 <f(t-Tx 2 )dT. 

Since <p is bounded, by the Dominated Convergence Theorem, we have that 

e-±T—'dT = V {t), 

according to (2.5). This proves the continuity up to the boundary of the solution U, con¬ 
cluding the proof that the function U defined in (2.8) is a continuous solution to the problem 
(1.5). □ 


3. Relation with the Marchaud fractional derivative 


We prove here the relation between the parabolic equation studied in Subsection 2.2 
and the Marchaud fractional derivative. Namely, the Marchaud derivative is obtained, in a 
certain sense, as the trace operator of the extension given by the solution of (1.5). 


Proof of Theorem 1.1. By inserting the expression of U(x,t) from (2.9), we compute 


( 1 r°° 2 

—— / z 2 s e"^r~ s "V(t-T)dT-^(t) 

4 S 1 (s) J o 

Recalling property (2.1) of the kernel, we have that 


lim x 2s (U(x,t) — ip{t)) 

x->- 0 + 


lim — . 

x->o+ 4 s T(s) 


r°° „2 

/ ir 2 s e _TF T _s_1 

J o 


(<p(t - t) - ip{t)) dr 


lim — „ , . 
x->o+ 4 s T(s) 


f°° 

/ e -n- 

Jo tS+ 
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e~~ < 1 


Now 

and since p is bounded, we have that 

ilV-Azm < 2Mt— 1 e L‘ «l,oo)). 

On the other hand, recalling that p is locally C 7 we have that 

\<P(t) - <p(t - t)| < CT 7 . 

Hence, since 7 > s, 


Using the Dominated Converge Theorem, we obtain 

lim x~ 2s ( U (x, <) - p{t)) = / lim 

X—> 0 + 4 S 1 (s) J Q x—v 0 ‘ 


4 s r( s ) y 0 

Hence for c s = 4 s T(s), we obtain 

/“OO 

—c s lim x _2s {U{x,t) — p(t)) = / 

no+ J 0 

by definition (1.2). This concludes the proof. 


I iptf-ij-pif) 

p 4 r - 

-10+ T s+1 

- r) - y?(t) 

7-S+l 


pit) - - t) 


pS + 1 


= D V(t) 


(3.1) 


□ 


4. Applications: A Harnack inequality for Marchaud-stationary functions 

In this part of the paper we prove a Harnack inequality for functions that have a vanishing 
Marchaud derivative in a bounded interval J, namely we prove here Theorem 1.3. At this 
purpose, we use a known Harnack inequality for degenerate parabolic operators, that can 
be found in [ 6 ], see Theorem 2.1. For the reader convenience we recall in Proposition 4.2 
this result in the case n = 1 . 


4.1. Preliminary notions. We would like to point out that the result given in [ ] was 
introduced for n > 3. Nevertheless the same proof works also for n = 1 with some adjust¬ 
ments. We recall here the main hypotheses we need to apply the Theorem, adapted in our 
case n = 1. It is worth to say that this problem has been studied in a more general fashion 
in [ 11 ] and [ 12 ]. 

The degenerate parabolic 


du 

w(x) Si 


d_ 

dx 



(4.1) 


is given in Q = (— R, R) x (0, T), for R > 0. 

In this particular case, the conductivity coefficient (i.e. the coefficient in front of the x 
derivative) and the specific heat (the coefficient of the t derivative) coincide. In [ 6 ], the 
equation was studied for different coefficients. A more general form of the equation is given 


in these terms: 


du 

Mx)- 


d_ 

dx 



(4.2) 


In that case, one requires that 
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and an integrability condition (also known as a Muckehoupt, or A 2 weight condition) on the 
weight w, given by 


t(pt L w(x)dx ) (pi Lw> dx 


= c 0 < 00 , 


(4.3) 


for any interval J C {—R, R). The constant Co is indicated as the A 2 constant of w. In our 
case, of course, we are left with the condition (4.3). 

In addition we consider the functional space 

W:= ju € L 2 {0,T-H x {J,w)) s.t. G L 2 (0,T;L 2 (J,r;))| . 

We denote here by L 2 (J, w), the Banach space of measurable functions u with finite weighted 
norm 

f • \ 1/2 

/ |u| 2 u> dx ) < 00 , 






by H x (J,w) the completion of C°°(J) under the norm 


\ u \\l ,W,J = 


(■ u 2 + | d x u\ 2 )w dx) 


\ 1/2 


and by Hq(J,w ) the completion of Cq° under the norm 


\ u \\l ,W,J = 


\d x ii\ 2 w dx 


1/2 


The time dependent Sobolev space L 2 (0, T; Hq (J, w)) is defined as the set of all measurable 
functions u such that 




:fi 


(4.4) 


\u(x, t)\ 2 w(x) dx dt \ < 00 . 

! (0,T) x J J 

In this setting, we introduce the notion of weak solution of problem (4.1). 

Definition 4.1. We say thatu G L 2 (0, T; H l ( J, w)) is a weak solution of (4.1) in J x (0,T) 
if, for every rj G W, such that rj{0) = iy(T) = 0 we have that 

. . / du dr) drj\ 

'j(x) - u—- dxdt = 0. 

'JX(O.T) \dxdx at J 

We have the next proposition (see Theorem 2.1 in [6]). 

Proposition 4.2. Let u be a positive solution in (—R,R) x (0 ,T) of (4-1) and assume that 
condition (4.3) holds, with constant cq. Then there exists 7 = y(co) > 0 such that 

sup u < 7 inf u (4-5) 

x-(# .£) (to+^f-.to+p 2 ) x-(f ) 

holds for to G (0, T) and any p such that 0 < p < R/2 and [to — p 2 ,to + P 2 ] C (O,? 1 ). 


Remark 4.3. The reader can easily imagine the general situation in any dimension as 
explicated in Theorem 1.2 in [6], where the coefficient a{x) in (4.2) is a matrix and the 
domains are cylinders. We have stated the Harnack inequality in (0,T). Nevertheless with 
a change of coordinates in space and time, we can always say that the Harnack inequality 
holds in any subset of (Ri,R 2 ) x (ti,t 2 ), where R \, R 2 ,17, t 2 G R. 
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4.2. Reflection of the solution. We consider here that D s <p(f) = 0 in an interval J. By 
taking the reflection U of the solution of problem (1.5), we prove that U is a solution in a 
weak sense of (1.5) across x = 0. 

It is useful to introduce a weak version of the limit lim x l ~ 2s d x u(x,t). In this sense, we 

x—>0 

have: 


lim x 1 2s -^—{x,t) = 0 


Definition 4.4. We say that in a weak sense 

,dU, 

x^o dx 

if and only if, for any r) £ W, such that 77(0) = r](T ) = 0 we have that 

lim [ x 1 - 211 ^- r] dt = 0. 


x —^0 


dx 


(4.6) 


Lemma 4.5. Let U: K.x [0,oo) — > R be a solution of the problem (1.5) such that, in a weak 
sense lim x 1 ~ 2s d x U(x, t) = 0. Then the extension 

a:—>0+ 


U(x,t) := 


U(x,t), (x, t) e[0, +oo) x (0, T) 

U(—x,t), ( x,t ) g(— oo,0) x (0,T) 


is a weak solution of 


dCM 1 2s U(x,t)) = d_ f. ,!_ 2 s dU(x,t) 
dt dx V dx 

in ( —R , R) x (0, T). 

Proof. We claim that the extension U is a weak solution of (4.8), hence that 


=0 . 


J (-R,R)x(0,T) 

We compute, integrating by parts 
rR 


dx dx dt 


fu 


*ldx)dt 
ox ox / 


/ R 1 2s ^-(R,t) rj(R,t) dt — lim f x 1 2s ^-rjdt 

J o x ~J o 


f 

/ 0 \-'o 


R 


.dU 


dx 


x* \ xj dx I dt 


dx 


R 


1-2 s 


dU_ 

dx 


{R , t)r](R , t) dt — 


io \J o 


i _ 2s dU , \ , 

x —— 7i dx dt, 

dt 


(4.7) 


(4.8) 


(4.9) 


where we have used the weak limit in (4.6) and the fact that U solves equation (4.8). In the 
same way, one obtains that 
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therefore, by summing up, 


l(-R,R)x(0,T) 


R 1 - 2s [^-(R,t)r](R,t) — ^-(—R,t)ri(—R,t) \ dt 


/ o 


dx 


dx 


1 1—2s 


-R 


dU , \ , 
— Vdx dt. 


Hence 


,- 2 . 9U^,_ 0 Sv\ dxdt 


( — R,R) X (0,T) 


dx dx 


dt 


R 1 2s i ^-(R,t)i]{R,t) - ^(-R,t)r}(-R,t )) dt 


io 


dx 


rR 


l-R 


dr] 


~KT d ~ U ~al \ dx] dt 


dt 


dt 


R 1 

f 

l-R 


2S dt 
^1 2s T)r](x, T) — U(x,0)r](x,0) S J dx 


= 0, 


since r](x,T) = r/(x,0) = 0 and r](R,t) = r/(—R,t) = 0. This is the claim in (2.5), and we 
conclude the proof of the Lemma. □ 


4.3. The Harnack inequality for the Marchaud derivative. We show here that the 
Harnack inequality for the Marchaud derivative can be deduced from the Harnack inequality 
applied to the extension operator. 


Proof of Theorem 1.3. We consider U to be the extension of <p, as introduced in Theorem 
1.1. Since <p is nonnegative, given the explicit solution U in (2.8), the function U is also 
positive. Now, we reflect U and obtain U, as we have done in Lemma 4.5. 

We take at first J = (0, T). Since D s (p(t) = 0 in (0,T), we have by definition that 


lim x 2s 

tC —>-0 


dU(x, t ) 
dx 


= 0, 


and thanks to Lemma 4.5, we obtain that U is a weak solution of (4.8) in, say, (—R,R) x 
(0,T) for a given R > 0. Moreover, the function Ixl 1 ^ 2 ® satisfies the condition (4.3), and 
according to Proposition 4.2, we have that 


sup u < 7 inf u 

(to-¥.*o-$)x-(-fi,$) (to+^,to + l)x-(f,§) 

for any p such that 0 < p < R/2 and [to — p 2 , to + p 2 ] C (0, T). It suffices now to slice the 
domain at x = 0 to obtain that 


to- 


sup 

3p r 


,t 0- 


zi(0, t) < 7 inf u(0 , t), 

(to + ^-.to + P 2 ) 
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hence 


sup ip(t) <7 inf (p(t) 

(to + ^,tc+P 2 ) 


for any p such that 0 < p < R/2 and [t 0 — p 2 , t 0 + p 2 ] C (0, T). 

Now, in order to prove that the Harnack inequality holds on any interval J C 1, one 
considers a translation of U , namely for any 9 G R, the function Ug(x, t) := U(x,t + Q), and 
reflects it as Lemma 4.5. Then Ug is a weak solution of (4.8), and Ug(0,t) = ip(t + 9). One 
obtains then, as a consequence of the Harnack inequality for the solution Ug, the following 

sup ip(t + 9) < 7 inf ip(t + 9) 

(to-2£,to-4) (*o+^4o+P 2 ) 

for any p such that 0 < p < R/2 and [to — p 2 , to + p 2 } C (0, T). Therefore 

sup ip(t) < 7 inf (p{t) 

(t 0 -e-2£,t 0 -8-4.) (to-e+^,to-e+P 2 ) 

for any p such that 0 < p < R/2 such that [to — 9 — p 2 ,to~ 9 + p 2 ] C (0,T). As 9 is arbitrary, 

one concludes that 

sup ip(t) < 7 inf tp(t) 

(*o-2£,*o- 4) (to + ^,to+P 2 ) 

for any p such that 0 < p < R/2 and [to — p 2 ,to + p 2 } C J. This concludes the proof of the 
Theorem. □ 

Remark 4.6. We would like to point out that the Harnack type inequality obtained in 
Theorem 1.3 can be equivalently stated as follows. Let us define for every 8 > 0 and for 
every t £ R the sets 

15 1 

I(t,8) =[t- ^<5,t + -<5], 


15 


7, 


1 {t,$) =[t- y S ’ t ~ 4^ 

=[t- + 

With this notation, the Harnack inequality affirms that for every I(r,8) C J 

sup <p < 7 inf ip. 

I+(t,S ) 


5. Backward equation 

We now consider the case of the right Marchaud fractional derivative, denoted by 
The following result is true: 

Theorem 5.1. Let s € (0,1) and ip: R —> R be a bounded, locally C^ function for s < 7 < 1. 
Let U- : [0,00) x R — > K. be a solution of the problem 

dU(x,t) 1 — 2 sdU(x,t) d 2 U{x,t) 


dt 

17 ( 0 , t) = <p(t), 
lim U(x, t) = 0 . 


dx 2 


x 


dx 


(x, t ) G (0,00) x R 
t G R 


(5.1) 
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Then U— defines the extension operator for ip, such that 

T> s _cp(t) = - lim . c a x~ 2s {U-{x,t) - ip(t)), (5.2) 

where 

c s =4 s r(s). 

We do not repeat all the computations, that are very similar to the case of the left 
Marchaud-derivative D s = D^_. 

We only point out that if f7_ is a solution of (5.1), then [7_(x, t) = U(x,—t), where U 
is the solution of the differential equation in (1.5). Thanks to Theorem 2.2 and keeping in 
mind (2.9), we get 


1 /*oo 2 

U-(x,t) = 4s]? e~^T~ s ~ l ip(t + T)dT. (5.3) 

Recalling the computations in (3.1) and the properties of the kernel \I/ S , see formula (2.4), 
we obtain that 


2s r o° 


lim x “ s (C/_(x, t) - p{t)) = lim / 

x— >0+ x— >0+ 4 s i (s) J 0 

= lim 


- 2s e 4t t s 1 (p(t + t) — ip(t)) dr 


x— 10+ 4 s T(s) J 0 


00 <p(t + t) - ip{t) 
e 4t — 


rS+1 


Thus, using the same argument as in the proof of Theorem 1.1, we conclude that 


lim x s ( U-(x,t ) — ip(t)) = 


1 


cp(t + r) - <p(t) 


x— >0+ 


4 s r(s) 


i-S+1 


(5.4) 


(5.5) 


that is 


Dl</?(t) = — c s lim x s (U-(x, t) — <p(t )). 

a->0+ 

dp 


It is worth to say that Di s Df_</?(f) = ——. Hence, using a different notation we can write 


that 


D +‘P(t)= 37 T, DLp(t)= (-37 V- 


dt 


dt 


6. Appendix 


In the Appendix, we provide some details on the Marchaud derivative. 

First of all, as stated in the Introduction, the Marchaud fractional operator is well 
defined for p bounded, locally C^ for 7 > s. Indeed, we have that: 


f 


p(t) - p(t - t) 


r s+l 


dr = 


pit) - p{t - t) 


7-S+l 


dr - 


p{t) - pit - t) 


r s ~\-1 


dr. 


= h + h ■ 


Since p is bounded, we have 

h < 2||^||loo( R ) 


+ i dr - C S>VP . 


( 6 . 1 ) 


Moreover, p is locally Holder, hence in (0,1) we may write 

I pH) - vH~ t ) I < ct 7 . 
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Therefore 

h<C f T^- S ~ l dT<C s ^, 

Jo 

recalling that 7 > s. 

There are, in literature, many other definitions of fractional derivatives. The interested 
reader can consult, for instance, [17] or [20] for further details. Here, we recall only the 
Riemann-Liouville fractional derivative, defined as 


v° ± m = 


±1 


/(*=Ft) 


dr 


for s £ C, 0 < 5is 
2.4.17 and 2.4.18 in [14]), given by 

D a ± m := 


T(1 — s ) dt J 0 t- 
< 1, see [19], Definition 1.16 and the Caputo derivative (see formulas 


±1 


f(tTr) 


dr. 


r(i - s) J 0 t s 

The definitions of Caputo or Riemann-Liouville are related to the Marchaud definition. 
Indeed, as one can see in formula (13.2) in the monograph [20], the Marchaud derivative 
is an extension of Riemann-Liouville’s, with weaker conditions on the function /. For 
sufficiently smooth / (say absolutely continuous, for instance), integrating by parts in the 
Riemann-Liouville definition, one can deduce the Marchaud notion (see also Theorem 1.17 
in [19]). 

As a further remark, the Marchaud derivative coincides with the notion of fractional 
derivative given by Griinwald and Letnikov, see [10], and Theorem 20.4 in [20] for the proof. 

We like also to remark that, just adapting the constant c s given in Theorem 1.1 by fixing 
c s = , in (1.7) we straightforwardly obtain the definition (1.1). The advantage of this 

choice is that —> <p as s —> 0 + and D±tp —> ip' as s —> 1~. Indeed, it is well known 

that the Marchaud derivative is not defined for s = 0 and s = 1 because in those cases the 
integral term in (1.1) (as in (1.2)) does not converge. However, one is able to pass to the 
limit by using the constant term from definition (1.1), that in this sense plays a fundamental 
role. 
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